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ABSTRACT 

In the light of the question whether most early-type dwarf (dE) galaxies in clusters formed through 
infall and transformation of late-type progenitors, we search for an imprint of such an infall history in 
the oldest, most centrally concentrated dE subclass of the Virgo cluster: the nucleated dEs that show 
no signatures of disks or central residual star formation. We select dEs in a (projected) region around 
the central elliptical galaxies, and subdivide them by their line-of-sight velocity into fast-moving and 
slow-moving ones. These subsamples turn out to have significantly different shapes: while the fast 
dEs are relatively fiat objects, the slow dEs are nearly round. Likewise, when subdividing the central 
dEs by their projected axial ratio into flat and round ones, their distributions of line-of-sight velocities 
differ significantly: the flat dEs have a broad, possibly two-peaked distribution, whereas the round 
dEs show a narrow single peak. We conclude that the round dEs probably are on circularized orbits, 
while the flat dEs are still on more eccentric or radial orbits typical for an infalling population. 
In this picture, the round dEs would have resided in the cluster already for a long time, or would 
even be a cluster-born species, explaining their nearly circular orbits. They would thus be the first 
generation of Virgo cluster dEs. Their shape could be caused by dynamical heating through repeated 
tidal interactions. Further investigations through stellar population measurements and studies of 
simulated galaxy clusters would be desirable to obtain definite conclusions on their origin. 

Subject headings: galaxies: dwarf — galaxies: elliptical and lenticular, cD — galaxies: structure 
— galaxies: evolution — galaxies: kinematics and dynamics — galaxies: clusters: 
individual: (Virgo) 



1. INTRODUCTION: THE VARIETY OF EARLY-TYPE 
DWARFS 

Early-type dwarf galaxies largely outnu mber all other 
galax y types in rich galaxy clusters (e.g., Binggeli et al. 
1985D, while they are less abundant on average in galax y 
groups, and almost absent in the field (|Gu et al.l l2006f) . 
Among dwarf galaxies, the fraction of those with early- 
type morphology is significa ntly larger in dynamicall y 
more evolved environments (jTrentham fe TulTvl l2002f ). 
Apart from this global trend, early-type dwarfs also 
show a pronounced relation with local galaxy d ensity 
within a given cluster or group: IBinggeli et alj (|19871 ) 
found that their number strongly increases w ith local 
densi ty, a nalogous to the correlation for giants (jDresslerl 
I1980T ). and lTullv fc Trenthaml (|2008l ) observed that early- 
type dwarfs strongly cluster around major E/SO galaxies. 
These findings have often been interpreted as indicat- 
ing that early-type dwarfs are actually formed through 
environmental processes, thereby transforming late-type 
into early-type objects. Such mechanisms, mainly gas 
stripping and tidal effects on structure and star for- 
mation, have been investigated in vario u s studies (e.g., 
iMoore et all Il998t Ivan Zee et all 12004 iSabatini et ail 
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12001 iBoselli et al.ll2008l ). with the main conclusion that 
one or more of these mechanisms must be responsible for 
the majority of cluster dEs. 

Despite the apparently simple structure of early-type 
dwarfs, a closer investigation of their characteristics re- 
vealed a surprising complexity. Disk signatures, s uch as 
bars and weak spiral arms (cf. Uerien etaTI 12000). were 
identifie d in a significant f raction of bright Virgo clus- 
ter dEs (|Lisker et al . 2006a, and refe rences therein) , con- 
firmin g earlier indications of disks (|Sandage fc Binggelil 
119841 : Binggeli fc Cameronl [l99l . The distribution of 
projected axial ratios reveals that these objects are not 
sph eroidal, but are rat her shaped like thick disk galax- 
ies (|Lisker et al.l l200"7f ) . Similar shapes, though some- 
what ttnckeT J _are_Jound_^for__the brighter nonnucleated 
dEs ([Ferguson fc Sa ndage 1989), as well as for those with 
blue central co lors from ongoing or very recent residua l 
star formation (jVigroux et al.lll984t [Lisker et al.ll2006bf ). 
These dE subclasses do not show strong clustering, but 
are distributed mo re like the late-type cluster galaxies 
(|Lisker et a l. 2007), lending support to the above forma- 
tion scenarios. 

Only the nucleated dEs that do not show disk features 
or blue cores appear to keep up the image of dEs as 
spheroidal objects consisting of old stars, concentrated 
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Fig. 1. — Selection process and resulting subsamples. Left: Positions of all Virgo cluster member galaxies (grey dots), of those 149 
dE(N)s for which heliocentric velocities are available (grey circles), and of selected fast-moving and slow-moving dE(N)s (grey triangles and 
black squares, respectively; see the middle panel). The positions of the major elliptical galaxies are indicated by crosses (M49, M84, M86, 
M87, counterclockwise from South). The central region (see text) is shown as grey-shaded area. Middle, top: Distribution of heliocentric 
velocities of all dE(N)s (light grey), central fast-moving dE(N)s (dark grey), and central slow-moving dE(N)s (black). The black vertical 
lines mark the semi-interquartile range (±320 km/s) of velocities relative to the median (1228 km/s), calculated from all dE(N)s. The 
grey vertical lines mark five semi-interquartile ranges, beyond which objects are excluded. Middle, bottom: Distribution of projected axial 
ratios for the fast and slow subsamples. Right, top: Distribution of projected axial ratios of central dE(N)s. We subdivide them at an axial 
ratio of 0.8 into a flat and a round subsample. Right, bottom: Resulting distribution of heliocentric velocities for the two subsamples. 



in regions of high local density ((Ferguson &; Sandaed 
1198a iLisker et al.1 120071) . They ar e the dE subclass 
with the oldest st ellar populations (jLisker et al.1 120081 : 
IPaudel et al.1 [2009) , and were concluded to be poten- 
tial primordi al objects in a study of the Coma and For- 
nax clusters (Rakos & Sch ombertl 120041 ) . Indeed, semi- 
analytic galaxy evolution models tied to N-body simu- 
lations of hierarchical structure formation show at least 
qualitative consistency with the observed scaling rela- 
tions and colors of ea r ly-type dwarf and g i ant g alaxies 
(|De Riicke et al.ll200l Ijanz fc Lisker! [20081 12009 ). indi- 
cating that a "cosmological" formation of dEs from 3- 
sigma density peaks on top of the cluster potential is 
plausible as well. In the study presented here, our goal 
is to assess whether even the centrally concentrated nu- 
cleated dEs of the Virgo cluster show some indication of 
an infall history, linking their formation process to envi- 
ronmental effects. 

2. SAMPLE SELECTION AND DATA 

As central cluster region, we select an area encompass- 
ing the massive elliptical galaxies M84, M86, and M87, 
requiring that the projected distance to M87, to M84, or 
to their connecting axis, is less than 1.27°, correspond- 
ing to 0.35 Mpc with m-M = 31 mag (Fig. [TJ left). 
This selection accounts for the fact that the dEs are not 
only clustered aroun d M87, which marks th e intraclus- 
ter gas mass center (jBoehringer et all [19 94) . but their 
distri bution peaks in the region between M84/M86 and 
M87 (|Binggeli et al.l [1981 . 

Bas ed on the study of dE subclasses by ILisker et al.l 
(|2007f), we select only nucleated dEs (dE(N)s) that 
display no disk fe atures or blue central colors (cf. 
ILisker et al.l l2006allbl ). Only certain cluster members 



(jBinggeli et all 119851 119931 ) brighter than m B < 18.0 
mag are taken into account. This is the same magni- 
tude limit up to wh ich the Virgo cluster catalog (VCC, 
IBinggeli et aLl[l985l) was found to be complete. Our clas- 
sification of nucleated and non-nucleated dEs relies on 
the VCC. However, many apparently non-nucleated dEs 
actually host a faint nucleus hardly detecta ble with the 
imag i ng data on w hich the VCC was based (|Grant et al.l 
[20051 : iCote et al.l [20061 ). Nevertheless, using the VCC 
classification yielded the said differences in structure, 
color, and distribution of the subclasses, which is why 
we stick to this definition (also see the discussion in 
ILisker et afll2007l) . 

Heliocentric velocities are available for 46 dE(N)s 
in the central region through the NASA/IPAC Extra- 
galactic Database (NED), of which 29 are provided 
by the Sloan Digital Sky Survey (SP SS) data re- 
lease 4 (I Adelman- McCarthy et aLll2006f), and 17 from 
other sources (IBinggeli et alj|1985lll993t iConselice et al.l 
I2lii0ll ISlrnien & Prugniell 120021 : lEvstigneeva et al.ll2007l) . 
Outer axial ratios are measured by fitting an ellipse to 
the r-band isophote at a semimajor axis a = 2ai 1 i, r , 
wher e aya, r defines the half-light aperture (Lisker ct al. 
120071) . Ultraviolet-optical colors (Kim et al. 2009, 
in preparation) are obtained by combining total far- 
ultraviolet (FUV) and near-ultraviolet (NUV) magni- 
tudes provided by GALEX (Galaxy Evolution Explorer, 
iMartin et al.ll2005f) w i th tot al optical magnitudes, mea- 
sured by_^iske^Fayj2002)_on SDSS data release 5 im- 
ages (jAdelman-McCarthv et alj|2007l ). 

3. RESULT: GALAXY SHAPE VERSUS VELOCITY 

From the heliocentric velocities, we calculate relative 
line-of-sight velocities, using the median value rj = 1228 
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Fig. 2. — Map of axial ratio differences. Left: Map of the 
Virgo cluster, with the greyscale level indicating the difference be- 
tween the local median axial ratio of slow and fast galaxies, ac- 
cording to the scale given below the figure. In bright regions, slow 
galaxies are rounder than fast ones, while this is reversed in dark 
regions. The map was calculated using circular areas of 1.27° ra- 
dius (corresponding to 0.35 Mpc with m — M = 31 mag), sampled 
every 0.127° in right ascension and declination. This strong over- 
sampling leads to a smooth map. The velocity separation is done 
as in Fig. [T] (middle). A value is calculated only if there are at 
least three of each of slow and fast galaxies. Right: Comparison 
map indicating how many galaxies are included in the local circular 
areas for which the values in the left panel were calculated. 

km/s of all 2 Virgo cluster dE(N)s as zeropoint. We then 
use the semi- interquartile range £ = 320 km/s of all 
dE(N)s to divide the central dE(N)s into fast-moving and 
slow-moving ones with respect to the line of sight, accord- 
ing to their relative velocities (Fig. [IJ top middle) , not 
distinguishing between the direction of movement (pos- 
itive or negative velocities). Two galaxies with veloci- 
ties beyond 5£ are excluded, leaving 27 fast and 17 slow 
dE(N)s. 

The projected axial ratios of fast and slow galaxies turn 
out to be significantly different (Fig.[TJ bottom middle): 
while the slow ones have nearly round shapes, the fast 
ones are much more flattened. The median axial ratios 
are 0.90 and 0.64 for slow and fast galaxies, and the mean 
values are 0.84 and 0.63, respectively. A K-S test yields 
a probability of 0.1% that both subsamples are drawn 
from the same parent distribution. 

The existence of two distinct subgroups is confirmed 
when the selection process is reversed: we now divide 
the central dE(N)s based on their projected axial ra- 
tio (Fig. [IJ top right), and compare the resulting dis- 
tributions of heliocentric velocities. Again, the two sub- 
samples are clearly different. The round galaxies have a 
narrow, centrally peaked velocity distribution, while the 
flatter galaxies have a much broader distribution, pos- 
sibly with a two-peak structure (Fig. [TJ bottom right). 
The velocity dispersion, calculated as standard deviation 
with one clipping at ±2.3u (±1% for a normal distribu- 
tion), is 737 km/s for the flatter subsample, and only 373 
km/s for the round subsample. A K-S test on the dis- 
tribution of velocities relative to the median, considering 
only absolute values, yields a probability of 0.1%. 



4. DISCUSSION 



4.1. 



Subpopulations and their distribution 

Following IConselice et al.l (|2001l ). the comparably high 
velocity dispersion of the flatter dE(N)s (737 km/s) could 

2 Velocities are available through NED for 149 of 210 dE(N)s 
that are certain cluster members. 



be seen as typical for a galaxy population with an infall 
history; spiral and irregular galaxies hav e similar disper- 
sions (776 and 727 km/s, respectively, IConselice et al.l 
120011) . However, there is one caveat to this straight- 
forward interpretation: our galaxies were selected to 
be in the (projected) central cluster region, and there, 
the velocity dispersion of all gal axy types was found 
to increase (IConselice et al.ll200lD . The conclusion of 
IConselice et al.l was that most galaxies are on radial or 
highly eccentric orbits. In the light of these observations, 
the low velocity dispersion of the round central dE(N)s 
is even more remarkable (37 3 km/s, as compared to the 
462 km/s of giant ellipticals. IConselice et al.l l2001). Due 
to the selection of galaxies in a projected central region, 
it is the most plausible explanation that those dE(N)s 
moving slow along the line of sight are on nearly circu- 
lar orbits, whereas most of those with high line-of-sight 
velocities should be on more eccentric orbits. 

If this was the case, we should see, to some extent, 
a reverse effect when selecting galaxies in the projected 
outskirts: those on more circular orbits should now, on 
average, be the faster-moving ones along the line of sight. 
Indeed, when defining "slow" and "fast" locally, the clear 
axial ratio difference seen in the central area does not 
only disappear towards the outer regions, but is some- 
what reversed (Fig. [3]), supporting our interpretation. 
The circularized orbits of the round, slow dE(N)s might 
thus indicate that they are the first generation of dE s 
in the Virgo cluster (see also lBiviano fe P oggianti 2009). 
Based on simulated cluster g alaxies ( Springel et al.ll2005D 
in the semi-analytic model of lDe Lucia fc Blaizotl ( 2007). 
just above the mass range of dwarfs, we find indications 
that for galaxies in the projected cluster center, slow- 
moving galaxies have typically fallen in earlier than faster 
ones. 

While those dE(N)s could represent the low-mass con- 
tinuation of elliptical galaxies in the framework of ACDM 
structure formation dJanz fc Lisk er 2008, 2009 ]), thu s be- 
ing "primordial" dEs (cf. iRakos fc Schombertll2004t ). one 
could also imagine that they were formed at early epochs 
out of other galaxies by transformation processes. On 
the one hand, the observations appear to be partially 
consistent with th e galaxy harassment scenario (e.g., 
iMoore et al.l [1996), in which an infalling disk galaxy 
experiences several close encounters with massive sub- 
structures and is eventually transfo rmed into a dE. In 
N-body simulations of harassment ( Mastrop ietro et al.l 
2005), the overall shape of harassed objects becomes 
rounder with each further encounter. The rounder galax- 
ies could thus have experienced infall and transformation 
at earlier epochs, having resided in the cluster for a longer 
time, leading to more tidal interactions. On the other 
hand, a possible inconsistency appears when considering 
that the strongest harassment effects occur when galax- 
ies are orbiting within t he cluster center, causing signifi- 
cant dynamical heating ( Mastropietro et al. 2005). This 
would imply that the rounder galaxies should be orbit- 
ing in the (threedimensional) central region, deep within 
the Virgo cluster potential well. However, from an anal- 
ysis of a galaxy cluster in t he Millennium-II simulation 
(|Bovla n-Kolch in et al.| [2009) with a mass comparable to 
Virgo, we find no significant difference between the clus- 
tercentric radii of galaxies moving slow and fast along 



4 



T. Lisker et al. 



' 6 
* ■ ** * I 


r '* ; : 

^•■••'^ 


■ 


■ 

■ 9^ 




!■■■!■■ 

A : 

^ •: 

- ^ ; 

i 


■ 

■ 
■ 


■ 

A 

* • v 

A V A 



0.2 0.4 0.6 0.8 1.0 -14 -16 -18 

Projected axial ratio M r 



Fig. 3. — Ultraviolet-optical colors versus axial ratio. Top 

left: Projected axia l ratio versus g — i color (measured within the 
half-light aperture. ILisker et al.l l2008T) for all dE(N)s (grey dots), 
the flat subsample (grey triangles), and the round subsample (black 
squares). The vertical dotted line denotes the subsample division 
at an axial ratio of 0.8. The circle and error bars denote the mean 
and standard deviation for each subsample. Top right: Color- 
magnitude relation. The solid grey line corr esponds to th e linea r 
fit for the high-density dE(N) subsample of ILisker et alj ( 2008). 
Bottom panels: Same as top panels, but for NUV — r color, using 
total magnitudes for color calculation. Only objects detected by 
GALEX in NUV are shown. For absolute magnitudes we adopt 
m — M = 3f mag. 

the line-of-sight. 3 Moreover, rounder galaxies end up 
with smaller sizes on average in the simulations, starting 
from the same progenitor mass. We do not observe such 
a difference in size; only the scatter is somewhat larger 
for the rounder/slower than for the flatter /faster ones. 

4.2. Stellar population aspects 

If the flatter dE(N)s experienced more recent infall, 
their star formation would have been truncated later, 
which should be reflected in their stellar population char- 
acteristics. Indeed, the g — i colors of the round dE(N)s 
are redder than those of the flat ones (Fig. [31 top left). 
This is confirmed when examining the color-magnitude 
relation (CMR, top right). The difference is even clearer 
when the age-sensitive NUV — r color is used instead 
(Fig- El bottom). A K-S test on the color residuals about 
the CMR yields a probability of 5.0% in g - i and 2.1% 
in NUV — r for the same parent distribution of the two 
subsamples. The values change to 7.4% and 3.7% when 
the reddest galaxy, which is an outlier to the CMR, is 
omitted. In FUV — r (not shown), where useful mea- 
surements are only available for ~ 50% due to detection 
limits, the colors of both subsamples show no difference 
within the errors. 

This behavior would be consistent with a small 
(~1 Gyr) age difference between two po pulations that 
are alr eady relatively old: as illustrated in ILisker fe Hani 
(2008), evolutionary synthesis models show that an ag- 
ing stellar population continues to become redder with 

3 Unfortunately, good -quality individu al distances from surface 
brightness fluctuations (Mci et al. 2007) are only available for 3 
slow and 2 fast galaxies. 



time in NUV — r, while it reaches a more or less constant 
FUV — r value. In fact, when fitting SDSS ugriz magni- 
tudes, pa rtly combined with near-infrared H band pho- 
tometry (|Gavazzi et al.l 120031 with improved reduction 
by T. Lisker), to GAL EV evolutionary synthesis mod- 
els l|Kotulla et all [2009). more than 60% of round/slow 
dE(N)s are fit best with an exponentially declining model 
with t = 0.5 Gyr. 4 In contrast, the best-fit model of 
more than 60% of flat/fast dE(N)s has a constant star 
formation rate for > 2 Gyr, followed b y an exponential 
decline with r > 1 Gyr (Hi elscheril2009l ). 5 

So far, we neglected the fact that the relation between 
axial ratio and color is not only seen for our sample 
of central dE(N)s, but is actually present, to some ex- 
tent, for the whole population of V irgo dE(N)s (Fig. [31 
grey d ots). This is not a new result: iBarazza fc Binggelil 
(2002) reported this very correlation, and concluded that 
it is a metallicity effect. Since outflows of metal-enriched 
gas occur primarily along the minor axis, flatter galax- 
ies have a lower metallicity. Nevertheless, the star- 
formation history and starburst episodes determine the 
gas and metallicity los s , respectively, from dwarf galax- 
ies (jRecchi et al.l 120061: IValcke et al] I2008Q . It therefore 
remains unclear whether these two findings, the corre- 
lation between shape and velocity on the one hand and 
between shape and color on the other hand, are at all 
related, or whether the latter correlation is simply inher- 
ent to all dEs, regardless of their velocity characteristics. 
More high-quality spectroscopic data are needed to in- 
vestigate whether differences in age, metallicity, and/or 
other stellar population characteristics are the primary 
cause of the observed color differences. 

As an aside, we note that a relation might exist be- 
tween the nucleus "strength" (i.e., its relative light con- 
tribution) and the host galaxy's projected axial ratio: the 
strongest nuclei were found by Binggeli et all (j2000D to 
res ide in nearly round dE(N)s. However, a similar study 
bv lGrant et al.l (|2005l ). now using Sersic instead of King 
profiles for the host galaxies, did not confirm a significant 
correlation except for a mild tendency 

5. CONCLUSIONS 

We have established that different populations, or gen- 
erations, of Virgo cluster dE galaxies can be distin- 
guished based on their projected shapes and line-of-sight 
velocities. This is most likely interpreted with the flat- 
ter objects moving on radial or highly eccentric orbits, 
indicating an infalling galaxy population. In contrast, 
the round dEs have circularized orbits, characteristic for 
cluster-born galaxies or very early infall. They might 
thus be the first generation of Virgo cluster dEs. While 
more detailed stellar population measurements would be 

4 The model grid comprised exponentially declining models with 
various decay times r (0.5, 1, 2, and 4 Gyr), as well as models with 
constant star formation rate until an age 8 of 2 or 8 Gyr, followed 
by an exponential decline with r of 1 or 3 Gyr. Models with fixed 
metallicity ([Fe/H] of 0.0, -0.3, -0.7 or -f .7) as well as chemically 
consistent models were calculated. The fits were performed by 
the dedicated x 2 minimization routine GAZELLE written by R. 
Kotulla (in preparation). 

5 A significantly different mean age of the two subsamples is, 
however, not found. Useful spectroscopic data are only available 
for 1 slow and 3 fast dE(N)s; a Lick index analysis yields an age of 
3 Gyr fo r one fast dE(N), a nd ages above 10 Gyr for the remaining 
objects (Paudcl et al. 2009). 
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necessary for an unambiguous interpretation, the obser- 
vational findings themselves constitute an important step 
towards understanding this most abundant galaxy popu- 
lation in clusters. Our results demonstrate the necessity 
to perform a similar deconvolution of dE populations in 
any ongoing and future study of dEs, if we aim at tracing 
back their evolutionary history along with that of their 
host clusters. 
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Fig. A.l. — SDSS images of fast and slow dE(N)s. Illustration of the selection by velocity shown in Fig.[T] The images were created 
by co-adding the SDSS g, r, and i bands. Contrast and scale is the same for all images. The angular scale is 1' X 1', corresponding to 
4.6 X 4.6 kpc 2 at m — M = 31 mag. 



